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ABSTRACT 


-V 

A  shipboard  cm  coherent  dual-polarized  focused  radar 
operating  at  9.23  GHz  and  an  optical  device  which  detects 
the  occurrence  of  specular  reflection  at  the  radar  incidence 
angle  have  been  used  to  investigate  mechanisms  governing 
radar  backscatter  froa  surface  water  waves  as  part  of  JOWIP 
(The  Georgia  Strait  Experiaent).  The  relative  contributions 
of  Bragg  and  specular  backscatter  aechanisas  have  been 
identified  using  the  output  of  the  optical  specular  detector 
and  the  radar  backscatter  polarization  ratios.  Analysis  of 
the  radar  data  together  with  data  froa  the  OREP  laser  wave 
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slope  gauge  Mhighes  and  Dawson. -this  issue)  has  been  used  to 
further  identify  the  aechanisas  governing  radar  backscatter 
and  aodulations  of  radar  backscatter  froa  surface  water 
waves.  It  is  found  that  specular  reflection  is  very 
significant  at  20^  incidence  angle  and  not  negligible  at 
40°.  It  is  found  that  when  specular  contributions  are 
absent  the  relationship  between  the  surface  wave  slope  and 
the  radar  backscatter  is  well-described  by  Bragg  theory. 
Specular  contributions,  when  they  occur,  appear  as  large 
spikes  on  top  of  the  Bragg  return.  Analysis  of  the  wave 
slope  data  shows  that  the  wave  slope  aodulations  at  "X-band" 
and  "L-band”  Bragg  wave  frequencies  are  comparable,  in 
contradiction  to  standard  wave  Modulation  theories  which 
predict  that  "X-band"  waves  should  exhibit  such  weaker 
aodulations.  Modulations  of  the  radar  backscatter  in  cala 
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seas  are  due  priaarily  to  Brags  scattering.  As  the  wind 
speed  increases,  the  hydrodynamic  modulations  tend  to 
decrease  but  the  radar  modulations  can  be  enhanced  by 
specular  reflection.  Possible  additional  contributions  due 
to  wedge-diffraction-type  scattering  at  higher  wind  speeds 
have  not  been  addressed  in  this  experiment.  Deep  minima, 
which  may  be  due  to  the  presence  of  organic  films  on  the 
surface,  have  been  observed  in  some  radar  returns  near 
current  patterns. 
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MICROWAVB  SCATTERING  FROM  INTBRNAL  WAVE  MODULATED  SURFACE 
WAVES:  A  SHIPBOARD  REAL  APERTURE  COHBRBNT  RADAR 

STUDY  IN  JOWIP  ( THB  GEORGIA  STRAIT  EXPERIMENT) 


1.  INTRODUCTION 

During  JOWIP  (The  Georgia  Strait  Bxperiaent),  synthetic 
aperture  radar  (SAR)  iaaging  of  internal  waves  was 
investigated  by  obtaining  SAR  iaages  of  internal  waves  at  X- 
and  L-bands  while  a  variety  of  aurface  and  subsurface 
aeasureaents  were  being  aade  to  provide  quantitative  data  on 
each  of  the  aechanisas  linking  internal  wave  displaceaents 
to  SAR  inage  intensity  variations  (see  Hughes  and  Dawaon, 
this  issue).  In  particular,  aeasureaents  were  aade  of 
surface  currents,  surface  wave  slopes,  real  aperture  X-band 
backscatter  froa  a  shipboard  radar,  the  occurrence  of 
specular  facets  (at  the  radar  incidence  angles)  on  the  water 
surface,  and  neteorological  data  such  as  wind  speed  and 
direction.  TRW  provided  the  shipboard  X-band  radar  and  a 
video  systea  for  detection  of  specular  facets.  This  paper 
describes  the  equipaent,  the  aeasureaents,  and  how  the  data 
has  been  used  together  with  the  wave  slope  data  to  identify 
the  aechanisas  responsible  for  radar  backscatter  and 
backscatter  aodulations. 

The  shipboard  X-band  radar  was  deployed  with  the 
following  two  objectives: 

(i)  To  deteraine  the  iaportance  of  specular  reflection  at 

20°  and  40°  incidence  angles  by  using  a  video  caaera  to 


detect  the  presence  of  specular  facets  while  aicrowave 
scattering  is  taking  place. 

(ii)  To  test  scattering  Models  by  aeaauring  the  in-situ  X- 
band  backscattering  cross-section  of  the  ocean  and 
comparing  the  cross-section  Modulation  with  the  output 
of  the  laser  slope  gauge  (deployed  by  B.  Hughes  of  the 
Defense  Research  Bstablishaent  of  the  Pacific,  or 
DRBP) . 

Neasureaents  to  address  the  above  two  objectives  were 
■ade  for  both  aabient  wind  wave  conditions  and  wind  waves 
which  were  undergoing  Modulations  by  internal  waves. 

2.  EXPBRIMBNT  DESCRIPTION 

Equipaent .  The  radar  except  for  ainor  Modifications, 
is  the  saae  unit  that  we  have  used  in  previous  laboratory 
wave  tank-wind  tunnel  acattering  studies . ( 1  -  2 >  It  is  an  X- 
band  CW  dual-polarized  systea,  operating  at  9.13  GHz  with 
roughly  1  watt  output.  The  antenna  is  a  22.9  cm  aperture 
conical  corrugated  horn  with  a  Matched  dielectric  lens  for 
spot  focusing.  The  3  dB  (2  way)  spot  size  is  18  ca  at  the  1 
a  focal  plane.  The  radar  is  a  superheterodyne  systea,  with 
the  transaitter  frequency  offset  froa  the  local  oscillator 
frequency  by  30  MHz.  Bach  polarization  channel  (vv  or  hh) 
has  two  aodes  of  output.  The  "raw  data"  aode  is  siaply  the 
30  MHz  IF  signal  downconverted  to  20  KHz.  This  signal  can 
be  used  for  accurate  power  aeasureaent  and  Doppler  spectrua 
analysis.  The  "processed  data"  aode  has  separate  linear 


amplitude  (i.e.,  not  power)  and  phase  (-180°  to  +180° ) 
output  channels  for  chart  display. 

For  detecting  specular  points,  a  bright  diffuaed  light, 
consisting  of  an  aircraft  landing  light  with  a  diffuser  in 
front,  is  aounted  on  one  side  of  the  horn  antenna  and  a 
video  caaera  is  aounted  on  the  other  side.  Both  are  pointed 
at  the  focal  spot  of  the  horn  on  the  water  surface.  A 
drawing  showing  the  aounting  arrangement  is  shown  in 
Figure  1. 

The  horn,  leap  and  video  caaera  were  aounted  onto  the 
frame  of  DRBP's  laser  slope  gauge  onboard  the  CFAV 
Endeavour,  with  a  distance  of  roughly  1  a  froa  the  water 
surface.  The  incidence  angle  of  20°  or  40°  was  set  at  the 
beginning  of  the  day.  The  horn  and  video  equipment  were 
always  pointed  along  the  axis  of  the  ship. 

On  the  bow  of  the  ship  was  aounted  a  separate  video 
caaera  for  recording  the  general  features  of  the  wave  field 
the  ship  was  going  through. 

Data  froa  the  radar  unit,  the  tiae  code  generator,  as 
well  as  froa  soae  DRRP  equipaent  (such  as  the  laser  slope 
gauge,  current  meter,  wind  aeter,  wave  height  gauge)  were 
recorded  on  two  analog  recorders.  The  two  video  cameras’ 
outputs  were  recorded  on  two  video  recorders. 

Data  collected  consists  of  119  runs,  each  run  being  a 
aore  or  less  continuous  recording  ranging  froa  a  minute  to 
aore  than  1/2  hour.  For  natural  internal  waves,  there  are 


66  runs  of  which  36  have  SAR  overpasses.  Roughly  82k  of  the 
runs  have  good  data. 

Procedures .  Bach  aorning,  calibration  and  alignment 
were  carried  out  on  the  deck  of  the  Endeavour  before  the 
equipaent  was  lowered  close  to  the  water  surface.  Alignaent 
of  the  leap  and  video  caaera  with  the  focal  spot  of  the  horn 
was  first  checked  with  the  help  of  a  teaplate  aounted  at  the 
horn’s  focal  plane.  The  incidence  angle  for  the  day  was 
then  set  relative  to  DREP's  slope  gauge  axis.  A  aetallic 
sphere  was  then  swung  through  the  focal  plane  of  the  horn  to 
calibrate  the  radar  unit.  After  the  equipaent  was  lowered 
into  position,  the  height  above  the  water  surface  was 
visually  checked  froa  the  bow  and  adjusted  accordingly. 

The  only  experiaental  difficulty  we  have  encountered 
was  the  alaost  constant  tuning  that  was  required  for  the 
radar.  In  a  laboratory  environaent,  once  the  radar  unit  was 
tuned  to  null  away  stray  reflections,  it  would  stay  tuned 
for  hours.  In  the  ocean,  the  horn  antenna  was  aounted  onto 
the  ORBP  slope  gauge  which  swayed  slightly  froa  side  to  side 
when  going  through  a  wave  field.  Since  the  horn  was 
connected  via  two  28'  waveguides  to  the  radar  unit,  the 
awaying  aotion  induced  phase  changes  in  the  return  signal, 
effectively  detuning  it  continuously.  In  a  rough  sea  at  a 
20°  incidence  angle,  the  backscattered  signal  is  auch 
stronger  than  this  detuning  noise.  In  a  cala  sea  at  a  40° 
incidence  angle,  the  aignal/detuning  noise  ratio  could 
becoae  quite  severe.  Fortunately,  our  receiver  has  an  IF 


output  which  records  the  raw  signal  down  converted  to 
20  kHz.  When  this  IF  signal  is  down  converted  to  40  Hz, 
detuning  appears  as  a  sharp  peak  at  40  Hz,  distinct  fron  the 
real  backscattered  signal,  which  is  usually  Doppler  ahifted 
by  ten's  of  Hz  fro ■  40  Hz. 

3.  RBSULTS  OF  SPBCULAR  REFLECTION  MEASURBMBNTS 

Specular  reflection  occurs  when  a  point  on  the  water 
surface  has  its  surface  nornal  pointing  directly  in  the 
■icrowave  incidence  direction.  When  this  occurs,  back- 
scattered  power  increases  drastically  for  a  short  duration 
ao  that  the  output  is  like  a  spike.  By  its  nature,  specular 
reflection  has  no  polarization  dependence  (in  contrast  to 
Bragg  scattering),  so  that  vv  and  hh  scattering  have  equal 
power*.  These  two  characteristics  can  be  used  to  identify 
specular  reflection  events  in  aicrowave  backscattering.  For 
positive  identification  of  specular  reflection,  however,  we 
have  to  rely  upon  the  coincidence  of  bright  patches  or  dots 
on  the  video  aonitor  with  spikes  in  the  aicrowave  beck- 
scattered  power.  Figure  2  shows  an  exaaple  of  this.  The 
top  row  shows  three  different  pictures  on  the  TV  aonitor  at 
three  different  tiaes.  It  can  be  seen  that  the  presence  of 
bright  dots  or  patches  coincides  with  occurrence  of  high 
outputs  (spikes)  in  both  the  vv  and  hh  channel. 

Furtheraore,  vv  and  hh  have  the  saae  power  level  in  the 
spikes.  Between  the  two  specular  events  there  are 

*Where  vv  and  hh  refer  to  the  polarization  states:  vv  is 
vertical  transait-vertical  receive  and  hh  is  horizontal 
transait-horizontal  receive. 


nonapecular  event*  where  vv  aaplitude  i*  roughly  twice  that 
of  hh,  aa  expected  froa  Bragg  acattering  at  40°  incidence 
angle. 

Trace*  3  and  4  ahow  DHEP’a  alope  and  wave  height 

Bbsaureaent  for  the  aaae  tiae  period.  The  occurrence  of 

apecular  event*  cannot  be  inferred  in  any  way  froa  theae 

trace*.  One  reaaon  for  thia,  of  courae,  ia  that  alope  and 

wave  height  gauge*  are  not  co-located  but  are  about  0.7  a 

away  froa  the  radar  antenna  footprint.  The  aajor  reaaon, 

however,  ia  that  the  alope  and  wave  height  aeaaureaenta  are 

point  aeaaureaenta  whereaa  the  radar  backacattering  ia  an 

area  aeaaureaent.  For  detecting  an  event  which  ia  aaall  in 

apatial  extent  and  ahort  in  tiae,  auch  aa  the  apecular 

facet,  an  area  aeaaureaent  will  ahow  a  auch  higher  frequency 

of  detection  than  a  point  aeaaureaent.  Suppoae  r«  =  radiua 

of  radar  antenna  footprint,  rap  -  radiua  of  curvature  of 

apecular  point,  r£  =  radiua  of  laaer  beaa  on  the  water 

aurface.  If  r«  >>  r*p  >>  r^ ,  it  i*  eaay  to  ahow  that 

probability  of  detecting  a  apecular  point  by  laaer  beaa  (hs\z 
probability  of  detecting  a  apecular  event  by  radar  \r2/ 

Since  the  10  dB  radiua  of  the  antenna  =  15  ca,  the  apecular 

point  radiua  of  curvature  typically  st  a  few  aa  and  the  laaer 

beaa  radiua  0  ( 1  aa),  the  above  probability  ratio 

2 

»  0(iy)  -  0(900). 

Thia  ia  the  reaaon  why  the  apecular  reflection  event* 
detected  by  the  radar  can  only  be  correlated  with  the  video 
recording  (an  area  aeaaureaent)  rather  than  the  alope  gauge 
or  wave-height  gauge. 
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To  deteraine  the  aaount  of  specular  reflection  under 
different  aea  states  and  incidence  angles,  we  use  the 
following  setup.  The  video  tape  containing  the  video  caaera 
recording  is  replayed  onto  a  TV  aonitor.  A  caaera  with  a 
photodiode  in  its  fila  plane  views  the  TV  aonitor  con¬ 
tinuously.  The  photodiode  output  will  show  spikes  whenever 
a  bright  patch  or  dot  appears  on  the  TV  screen.  The  spikes 
occur  in  groups  and  within  the  groups  they  are  roughly  1/60 
sec  apart,  1/60  sec  being  the  TV  field  rate.  The  aagnitude 
of  each  spike  is  proportional  to  the  aagnitude  of  the 
specular  return.  The  analog  tape  containing  the  aicrowave 
data  is  also  replayed  siaultaneously .  If  the  video  tape  and 
the  analog  tape  are  synched  correctly,  the  photodiode  output 
spikes  should  coincide  with  a  specular  event  in  the 
aicrowave  channel.  Since  both  video  recorder  and  analog 
recorder  are  aechanical  devices,  synching  both  to  the 
accuracy  involved  (1/10  sec)  over  a  long  period  of  tiae  is 
alaost  impossible.  We  solved  this  problea  by  inputting  both 
photodiode  output  and  the  aicrowave  data  into  a  digital 
scope  with  large  aeaory.  Proper  synching  is  then  perforaed 
on  the  digital  scope.  Using  this  aethod,  roughly  10  ainutes 
of  data  can  be  analyzed  at  a  stretch.  For  each  continuous 
segaent  of  data,  the  total  aaount  of  power  due  to  specular 
reflection  as  well  as  the  total  tiae  during  which  specular 
reflection  occurs  are  suaaed.  (We  are  assuaing  for 
siaplicity  that  whenever  the  video  caaera  sees  bright  dots 
and  patches,  the  corresponding  aicrowave  power  coaes  from 


specular  reflection.)  These  are  then  expressed  as  a 
fraction  of  the  total  power  and  a  fraction  of  the  total  tiae 
of  the  data  segaent.  The  results  as  a  function  of  wind 
speed  and  incidence  angle  are  shown  in  Figures  3  and  4.  In 
order  to  obtain  as  statistically  aeaningful  a  set  of  results 
as  possible,  about  52  ainutes  of  data  were  processed  froa  a 
total  of  21  runs  covering  a  range  of  wind  speeds  froa  0.7  to 
8.6  a/sec  and  two  incidence  angles  (20°  and  40°).  Except 
for  the  cases  as  indicated,  all  the  data  points  are  froa 
pure  wind  wave  runs,  i.e.,  where  internal  waves  are  absent. 
In  all  cases,  the  ship  is  sailing  directly  into  the  wind 
130°.  Froa  the  figures,  we  can  draw  the  following 
conclusions : 

(i)  For  20°  incidence,  specular  power  as  a  fraction  of 
total  power  varies  froa  —  70*  at  low  wind  (~  1  a/s)  to 
90%  at  high  wind  (~  8  a/s). 

(ii)  For  40°  incidence,  specular  power  increases  froa  alaost 
nothing  at  low  wind  to  ~  30-40%  at  7  a/s. 

(iii)  For  20°  incidence,  the  tiae  during  which  specular 
reflection  occurs  is  20%  of  the  total  tiae. 

(iv)  At  40°  incidence,  the  presence  of  an  internal  wave 
field  increases  power  in  specular  reflection  quite 
noticeably.  This  probably  is  a  result  of  enhanced 
induced  "wavebreaking"  by  the  internal  wave. 

We  want  to  point  out  that  the  above  conclusions  are 
froa  data  sets  obtained  in  Oabob  Bay  and  Boundary  Pass  where 
the  fetch  is  highly  liaited.  It  does  show,  however,  that 


specular  reflection  at  X-band  ia  sore  significant  than  has 
been  generally  expected.  It  has  generally  been  aaauaed  that 
specular  contributions  should  be  negligible  at  incidence 
angles  greater  than  20°  (see,  e.g.,  Valenzuela,  1978). 

An  often  heard  consent  discounting  the  possible 
inportance  of  specular  reflection  at  X-band  ia  that  "facets 
on  the  water  surface  which  are  optically  specular  nay  not  be 
specular  to  nicrowave  radiation".  The  exact  aeaning  of  this 
stateaent  is  that  the  radius  of  curvature  at  the  specular 
points  nay  be  nuch  snaller  than  the  nicrowave  wavelength 
(but  would  still  be  nuch  longer  than  optical  wavelength). 

In  this  case,  nicrowave  scattering  fron  the  specular  points 
would  be  Rayleigh  scattering,  i.e., 

°*icro«»*»  «  trr*(kr>4 

where  r  =  radius  of  curvature  at  specular  point 

k  =  2ir/X 

X  =  nicrowave  wavelength 

For  optical  wavelengths,  scsttering  would  be  in  the 
geonetric  optics  region  so  that  cross  section  would  be  given 
by 

^optical  t  irr* 

Thus,  if  r  <<  X,  (kr)4  will  be  <<  1  and  cBlcroHsvt  << 

^optical  . 

To  answer  the  question  whether  specular  reflection  at 
X-  and  L-bands  is  predoninant ly  Rayleigh  scattering,  we  need 
to  know  the  size  distribution  of  the  specular  facets  (or, 


■ore  accurately,  the  radiua  of  curvature  dietribution  of  the 
apecular  pointa).  This  can  be  estimated  by  analyzing  the 
video  images  we  have  recorded.  Since  our  illumination 
source  is  a  finite  size  diffused  laap,  the  apparent  size  of 
the  iaage  of  the  laap  in  water  at  each  specular  point  will 
be  directly  related  to  the  radiua  of  curvature  at  that 
specular  point  (see  Appendix  A).  This  analysis  has  been 
carried  out  for  two  wind  wave  runs  at  40°  incidence.  Figure 
5  shows  the  size  distribution  of  the  facets  for  one  run 
plotted  as  histograms  in  kr  for  both  X-  and  L-bands. 
Generally  speaking,  kr  <  1  can  be  regarded  as  the  Rayleigh 
scattering  region,  1  <  kr  <  0(10)  can  be  regarded  as  the 
resonance  region  and  kr  >  0(10)  can  be  regarded  as  the 
geoaetric  optics  region.  Blooaing  on  the  TV  screen  aay  be  a 
significant  source  of  error  in  aeasuring  the  size  of  saall 
facets . 

Figure  6  shows  power  distribution  of  the  facets  plotted 
as  histograas  in  kr,  i.e.,  the  vertical  axis  shows  total 
power  contributed  by  facets  of  each  particular  kr  value.  It 
can  be  seen  that  for  X-band,  aost  of  the  specular  reflection 
occurs  in  the  resonance  and  geoaetric  optics  region,  i.e., 
the  facets  are  not  aaall  relative  to  the  X-band  wavelength 
at  40°  incidence.  For  L-band,  there  is  aore  Rayleigh 
acatterlng  but  still  aost  of  the  specularly  reflected  power 
occurs  in  the  resonance  region  and  geoaetric  optics  region. 
Our  aeasureaents  aay  be  subject  to  errors  due  to  blooaing 
(as  mentioned  above)  and  aotion  blurring.  Blooaing  aay 


iapose  a  lower  liait  on  size  aeasureaent.  Its  error  for  big 
patches  should  be  such  saaller  in  proportion.  He  have 
calibrations  giving  us  a  rough  estiaate  of  the  error  due  to 
blooaing.  For  aotion  blurring,  the  correct  solution  (which 
aay  be  iapleaented  in  future  tests)  should  be  the  use  of 
either  a  fast  shutter  in  front  of  the  video  caaera  or  a 
synchronized  flash  laap.  For  our  recorded  data,  we  can 
aeasure  only  the  cross-track  size  of  the  patches.  This 
presumably  is  not  subject  to  aotion  blurring.  With  all 
these  caveats,  we  think  it  is  still  safe  to  conclude  that  at 
X-band  at  40°  incidence,  specular  reflection  occurs 
principally  in  the  resonance  and  geoaetric  optica  region. 


The  relative  contributions  of  specular  and  Bragg 
backscatter  to  the  total  return  as  a  function  of  the  size  of 
the  radar  illumination  spot  on  the  water  surface  has  not 
been  investigated  in  these  experiaents.  A  question  can  be 
raised  as  to  how  the  relative  contributions  sight  vary  as 
platfora  altitude  and/or  spot  size  is  varied.  Specular 
contributions  would  add  incoherently  and  sight  be 
overwhelmed  by  Bragg  contributions  which  would  increase 
faster  than  linearly  if  they  continued  to  add  coherently. 
Since  water  waves  at  centiaeter  scales  are  expected  to  have 
correlation  scales  of  no  aore  than  a  few  of  their  own 
wavelengths,  we  expect  that  as  the  illumination  spot  is 
increased  over  that  used  in  this  experiment,  the 
contributions  frop  both  Bragg  and  specular  aechanisas  will 
add  incoherently  and  the  relative  contributions  should  be 
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close  to  what  we  have  observed  here  even  for  platforms  at 
■uch  higher  altitudes.  We  intend  to  investigate  this  aspect 
of  the  problea  quantitatively  in  future  experiments.  We 
wish  to  emphasize  that  the  video  system  was  meant  from  the 
beginning  to  detect  specular  points.  It  was  only  after  the 
data  had  been  collected  that  we  realized  some  semi- 
quantitative  information  could  be  gathered  about  specular 
reflection  in  general  and  about  the  radii  of  curvature  of 
specular  points  in  particular.  The  data  and  analysis 
presented  here  must  be  regarded  in  this  light.  A  definitive 
study  of  specular  reflection  will  have  to  wait  for  improved 
instrumentation  and  techniques. 

4.  COMPARISON  OF  MICROWAVE  BACKSCATTERING  MODULATION 

WITH  SURFACE  SLOPE  MODULATION 

The  DREP  laser  slope  gauge  offers  us  a  unique 
opportunity  to  test  models  of  microwave  scattering  from  the 
ocean  surface.  Although  the  slope  gauge  is  not  co-located 
with  the  microwave  beam  so  that  a  completely  deterministic 
study  of  scattering  (such  as  was  done  in  the  laboratory. 
Refs.  1,  2)  is  not  possible,  its  close  proximity  (73  cm  to 
the  left  when  looking  forward)  means  that  modulations  in 
microwave  backscattering  and  in  wave  slope  can  be  compared 
to  provide  some  information  on  the  validity  of  various 
scattering  models. 

We  will  first  discuss  how  the  microwave  and  slope  data 
are  analyzed.  Microwave  data  consists  of  two  separste 
channels,  vv  and  hh,  down-converted  in  two  steps,  from 


9.13  GHz  to  30  MHz  to  20  KHz  before  it  is  recorded.  To  do 
detailed  analysis,  the  recorded  signals  are  down-converted 
again  to  40  Hz  by  nixing  with  a  signal  roughly  at  20  KHz 
froa  a  frequency  synthesizer.  (A  synthesizer  is  required 
since  the  recorded  signal  deviatea  slightly  froa  20  KHz  froa 
day  to  day).  The  signal  is  first  low  pass  filtered  at 
300  Hz  before  being  recorded  onto  a  digital  tape  recorder  at 
512  Hz  sanpling  rate.  The  low  pass  filtering  aininizes  high 
frequency  fold-over.  The  digital  tape  is  then  analyzed  by 
our  PRIME  750  ainicoaputer .  A  1024  point  Feat  Fourier 
Transfora  with  a  Hanning  weighting  is  perforaed  on  2  seconds 
of  data  every  0.5  aec.  The  spectrua  typically  shows  a  sharp 
peak  at  40  Hz  corresponding  to  the  detuned  signal  and  a 
broad  peak  Doppler  shifted  froa  40  Hz,  corresponding  to  the 
backscattered  signal.  For  caaes  where  the  Doppler  peak  is 
clearly  aeparated  froa  the  40  Hz  detuned  signal,  the 
spectrua  is  integrated  froa  50  Hz  to  200  Hz  to  give  the 
backscattered  power  as  a  function  of  tine  at  0.5  sec  inter¬ 
val.  For  caaea  where  the  broad  Doppler  peak  hovers  around 
40  Hz  (which  happens  when  the  ship  and  the  waves  are  noving 
at  the  aaae  velocity),  the  portion  of  the  sharp  detuning 
peak  above  the  broad  Doppler  peak  is  siaply  subtracted  away 
before  integration  to  yield  the  total  power  of  the  Doppler 
broadened  aignal.  This  process  of  integrating  the  spectrua 
to  obtain  the  backscattered  power  is  perforaed  for  both 
polarizations  to  produce  the  polarization  ratio  (hh/vv)  as  a 
function  of  tine  at  0.5  sec  interval.  The  aean  Doppler 


shift  as  well  as  the  width  of  the  Doppler  spectrua  (at  10  dB 
below  peak)  can  also  be  obtained  by  further  processing  for 
every  spectrua.  i.e.,  at  0.5  sec  interval.  To  suaaarize, 
the  radar  signal  is  analyzed  to  obtain  at  0.5  aec  interval 
the  following  data:  vv  power,  hh  power,  aean  Doppler  shift 
and  width  of  Doppler  spectrua  of  both  polarizations.  On 
the  PRIME  750,  a  typical  10-ainute  run  will  take  about 
10  ainutes  coaputer  tiae  to  analyze. 

The  laaer  slope  gauge  output  is  a  voltage,  which  after 
calibration,  becoaes  slope  as  a  function  of  tiae.  The  slope 
signal  analysis  procedure  is  siailar  in  soae  aspects  to  that 
of  the  radar.  It  is  low  pass  filtered  at  300  Hz  and  then 
recorded  onto  a  digital  tape  at  512  Hz  saapling  rate.  The 
tape  is  then  analyzed  on  the  PRIME  750  which  perforas  a  FFT 
on  2  sec.  of  data  every  0.5  sec.  The  resulting  spectrua  of 
every  2  seconds  of  data  is  the  slope  spectrua.  We  recall 
that  the  slope  power  spectral  density,  S(f),  is  related  to 
the  waveheight  power  spectral  density,  Y(f),  by: 

S(f)  =  k*Y(f) 

Let  us  assuae  that  an  X-band  Bragg  wave  with  wavelength 
and  X*  and  wavenuaber  ks  is  aoving  towards  the  ship  with  a 
relative  velocity  v.  The  laser  beaa  will  detect  this  Bragg 
wave  as  a  disturbance  with  a  "Bragg  frequency"  fs  where 

ft  =  v/X,  (1) 

If  Bragg  acattering  theory  is  correct,  then  the  aicrowave 
backscatterad  power,  P,  should  be  proportional  to  the 
waveheight  spectral  density  at  k>  in  k  -  space,  which  is 


ITJ-WM  JTU  KVKV1TU  WWWW  WWir.TJVXUWIS™ OTWWWJRWWm^WWVWWWWJV'JVV^i 


proportional  to  the  waveheight  apectral  denaity  at  f • ,  in 
frequency  apace,  i.e., 

P  *  Y(k»  )  «  Y(  fi  )  =  S(f») 

Thus  modulation  in  backacattered  power  ahould  be 
proportional  to  modulation  of  the  alope  apectrum  at  the 
'Bragg  frequency”  f*  =  v/X».  In  a  relative  dB  acale, 
modulation  in  microwave  backacattered  power  should  be 
identical  to  modulation  of  the  alope  apectrum  at  the  ”Bragg 
frequency”  within  a  constant  additive  factor. 

The  strategy  of  teating  Bragg  scattering  theory  is 
therefore  quite  simple  in  principle:  Plot  the  microwave 
backacattered  power  (vv  or  hh)  in  a  dB  acale  as  a  function 
of  tine  as  the  ship  passes  through  an  internal  wave  field; 
aelect  a  narrow  band  around  the  ”Bragg  frequency”  in  the 
alope  apectrum;  integrate  the  alope  power  in  this  narrowband 
to  obtain  the  ”X-band-Bragg  slope  power”  in  a  dB  scale  as  a 
function  of  time.  The  2  plots  should  be  identical  if  Bragg 
scattering  theory  ia  correct. 

To  implement  the  strategy  requires  knowing  where  the 
”Bragg  frequency”  is  and  how  wide  a  band  we  should  integrate 
around  it.  From  Bq.(l),  the  ”Bragg  frequency”  is  deter¬ 
mined  by  the  encounter  velocity  between  the  ship  and  the 
Bragg  wave.  We  can  determine  the  encounter  velocity  with 
increasing  degree  of  accuracy  and  complexity  by  the 
following  three  methods.  (1)  The  simplest  thing  to  do  is 
to  assume  the  ship  is  moving  with  a  constant  velocity  (which 
we  obtain  by  averaging  the  ship’s  bow  current  meter  reading 


over  a  long  period  of  time  before  the  ship  enters  the 
internal  wave  field).  We  also  assuie  the  Bragg  waves  are 
either  moving  toward  or  away  from  the  ship  at  their  natural 
dispersion  phase  speed.  The  encounter  velocity  is  simply 
the  sua  of  the  2  velocities.  Consequently,  slope  analysis 
is  performed  by  setting  a  bandpass  filter  at  a  fixed 
frequency.  (2)  The  ship,  however,  is  not  moving  at 
constant  velocity.  In  fact,  the  ship's  bow  current  meter 
clearly  shows  it  to  be  varying  when  the  ship  goes  through  an 
internal  wave  field.  The  ship's  varying  velocity  can  be 
taken  into  account  in  determining  the  encounter  frequency 
with  the  Bragg  wave.  This  corresponds  to  a  slope  analysis 
with  a  tracking  filter  that  tracks  with  ship's  varying 
velocity.  This  obviously  complicates  the  data  processing 
considerably.  (3)  Finally,  we  should  realise  that  the 
Bragg  wave  phase  speed  is  highly  uncertain.  The  Bragg  wave 
may  be  advected  by  a  wind  drift  layer  (typically  ~  3*  wind 
speed);  it  may  be  advected  by  orbital  speed  of  gravity 
waves;  or  it  may  actually  be  a  locked  component  at  the  crest 
of  a  short  gravity  wave.  At  this  point,  just  when  it  seems 
quite  hopeless  that  we  can  determine  the  encounter  velocity 
accurately,  we  should  realise  that  the  encounter  velocity 
has  already  been  measured  —  by  our  Doppler  radar!  In  fact, 
a  little  arithmetic  will  show  that  the  "Ship-Bragg  wave 
encounter  frequency"  is  identical  to  the  radar  measured 
Doppler  shift  frequency.  The  slope  analysis  should 


therefore  involve  a  tracking  filter  that  tracks  with  the 
radar  Measured  Doppler  frequency. 

We  note  here  that  this  analysis  of  the  wave  slope  data 
assuaes  that  the  surface  waves  are  noving  along  the 
direction  of  the  line  of  sight  of  the  radar,  which  is  also 
the  direction  of  the  ships  heading,  since  the  radar  responds 
only  to  Bragg  waves  aoving  in  that  direction.  In  a 
Multidirectional  situation  it  is  possible  that  waves  aoving 
at  an  angle  to  the  radar  line  of  sight  could  contribute  to 
slope  data  at  the  Bragg  wave  frequency  even  though  they  are 
not  waves  at  the  Bragg  wavelength  aoving  in  the  direction  of 
the  radar  line  of  sight.  The  slope  gauge  is  therefore 
capable  of  responding  to  waves  approaching  at  a  wider  range 
of  angles  than  the  radar  and  this  could  affect  coaparisons 
between  the  radar  backscatter  and  wave  slope  data  bandpassed 
at  the  "Bragg  frequency".  This  effect  is  considered 
negligible  for  the  results  reported  here  for  several 
reasons.  First,  the  internal  waves  were  long-crested,  the 
aeasureaents  were  Bade  while  the  ship  was  heading  into  the 
wave  fronts  (i.e.,  heading  perpendicular  to  the  long 
crests),  and  visual  and  video  observations  indicated  that 
the  short  waves  were  aligned  in  the  direction  of  the 
internal  wave  propagation,  so  that  the  configuration  was 
essentially  unidirectional.  Second,  since  the  aajor 
contributions  to  the  relative  speed  of  the  short  waves  and 
the  ship  is  froa  the  ship  speed,  it  can  rather  easily  be 
seen  that  waves  at  an  angle  to  the  ship  heading  (and  radar 


line  of  eight)  which  would  contribute  slope  data  at  the 
frequency  of  Bragg  waves  aoving  along  the  radar  line  of 
sight  are  waves  having  shorter  wavelengths  than  the  true 
Bragg  waves.  If  the  surface  wave  slope  spectrua  decreases 
with  decreasing  wavelength,  this  contributes  further  to  the 
essentially  unidirectional  nature  of  these  aeasureaents . 
Thirdly,  it  is  the  aodulations  of  wave  slope  and  radar 
backscatter  that  are  being  coapared  here  and  not  the 
absolute  values.  Begardless  of  the  directional  dependence 
of  the  slope  spectrua,  aodulations  of  surface  waves  induced 
by  internal  waves  are  strongest  in  the  surface  wave 
coaponents  aligned  with  the  internal  wave  direction,  so  that 
the  aodulations  of  interest  tend  to  be  concentrated  in  the 
unidirectional  wave  coaponents. 

The  question  of  what  bandwidth  should  be  used  in  the 
slope  analysis  is  quite  interesting.  Since  the  aicrowave 
Doppler  bandwidth  is  tiae  varying,  we  aay  be  teapted  to 
think  that  if  slope  power  aodulation  does  turn  out  to  be 
identical  to  aicrowave  power  aodulation,  the  bandwidth  used 
in  slope  analysis  should  also  be  tiae  varying.  This 
reasoning  is  false  on  closer  analysis.  We  can  understand 
this  aore  clearly  by  referring  to  Figure  7  which  shows  a 
schematic  plot  of  the  dispersion  relation.  In  the  ocean, 
due  to  the  reasons  aentioned  above,  the  dispersion  relation 
is  no  longer  a  sharp  line.  Rather,  it  is  a  band.  If  power 
is  plotted  vs.  u)  and  k  in  a  3-D  plot,  the  waveheight  power 
spectral  density  is  no  longer  a  sheet,  but  a  ridge.  If  the 
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Bragg  scattering  theory  is  correct,  the  radar  aeasureaent  is 
a  cut  across  this  ridge  at  a  constant  k,  with  a  finite  Ak 
due  to  the  finite  spot  size  of  the  antenna.  The  cut  at 
constant  k  across  the  ridge  is  in  fact  the  Doppler  spectrua 
of  the  radar  signal.  The  slope  spectral  analysis  (sore 
accurately,  the  waveheight  spectral  analysis)  is  a  cut 
across  this  ridge  at  a  constant  01,  with  a  finite  Aw  to  be 
deterained.  Froa  the  figure,  we  can  see  that  coaparing  the 
integral  along  the  constant  k-cut  (corresponding  to  aicro- 
wave  backscattered  power,  assuaing  Bragg  scattering  theory 
is  valid)  and  the  integral  along  the  constant  u -cut  (the 
Bragg  slope  power)  is  strictly  speaking  valid  only  when  the 
dispersion  relation  is  a  line.  When  the  dispersion  relation 
becoaes  a  band,  the  2  integrals  are  integrating  over  2 
different  regions  (shown  as  shaded  in  the  figure).  However, 
it  is  easy  to  show  that  the  2  regions  have  the  saae  area,  to 
first  approxiaation .  For  reasonably  behaved  power  spectral 
density  "ridges",  the  2  integrals  reaain  identical.  Froa 
Figure  7,  the  answer  to  our  original  question  should  now  be 
obvious.  The  Aw  should  be  chosen  to  correspond  to  the  Ak  of 


the  radar  footprint,  i.e.,  Aw  =  (dw/dk)Ak  ,  where  (dw/dk)  is 
the  slope  of  the  dispersion  relation. 

Barlier,  we  asked  whether  the  slope  analysis  bandwidth 
should  vary  when  the  width  of  the  aicrowave  Doppler  spectrua 
varies  with  tiae.  Froa  Figure  7,  it  is  clear  that  the 
Doppler  spectrua  changes  width  as  a  result  of  the  dispersion 
relation  band  changing  width,  e.g.,  froa  the  solid  lines  to 
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the  dashed  lines  in  Figure  7.  The  2  integrals  along  con¬ 
stant  k  and  constant  to  reaain  equal.  Hence,  Ato  should 
reaain  constant  because  Ak  of  radar  is  constant,  a  reasoning 
which  is  not  quite  obvious  without  referring  to  Figure  7. 

At o,  however,  does  not  always  reaain  constant,  tfhen  the 
ship  is  going  through  an  internal  wave  field,  its  velocity 
changes.  Consequently  the  dispersion  relation  aeasured  froa 
the  ship  changes  slope.  As  we  aentioned  above, 

Aw  a  (i^)Ak  where  ^  =  slope  of  dispersion  relation, 

function  of  the  ship’s  velocity 
Hence,  Au>  should  track  with  the  ship’s  velocity. 

To  suaaarize,  the  slope  analysis  is  carried  out  to 
coapute  the  "X-band"  Bragg  slope  power.  The  "Bragg 
frequency"  tracks  the  aean  Doppler  frequency  as  aeasured  by 
our  radar.  The  bandwidth  for  doing  the  slope  power 
integration  tracks  with  the  ship’s  velocity  as  aeasured  by 
the  bow  current  aeter. 

The  above  analysis  is  perforaed  also  for  the  "L-band" 
Bragg  alope  power.  In  this  case,  we  do  not  have  a  ship¬ 
board  L-band  radar  to  tell  ue  the  velocity  of  the  L-band 
Bragg  waves  relative  to  the  ship.  For  siaplicity  we  assuae 
that  the  relative  speed  of  the  ship  and  the  L-band  waves  is 
the  sane  as  that  for  the  X-band  waves.  This  is  equivalent 
to  assuaiog  that  the  X-  and  L-band  waves  are  phase  locked 
rather  than  dispersive,  but  we  note  here  that  the  results  in 
thse  cases  turn  out  to  be  insensitive  to  which  assuaption  is 
Bade  because  the  doainant  part  of  the  relative  velocity  is 
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the  ship's  own  velocity.  The  Au>  for  L-band  slope 
integration  is  taken  to  scale  with  wavelength,  i.e.,  about 
1/7  that  of  X-band. 

The  video  signal  has  been  analyzed  as  aentioned  pre¬ 
viously  on  a  f raae-by-f raae  basis  to  provide  a  seai- 
quantitative  analysis  of  specular  reflection.  The  video 
signal  is  a  free-running  indicator  of  the  occurrence  of 
specular  points.  The  signal  is  first  low-pass  filtered  at  5 
Hz  to  reaove  the  60  Hz  spikes  (the  TV  field  rate).  It  is 
then  recorded  at  512  Hz  saapling  rate  on  a  digital  recorder. 
The  power  is  integrated  over  a  2  second  interval  every  0.5 
asec.  The  logaritha  of  the  power  is  taken  before  it  is 
plotted  out  together  with  the  rest  of  the  signals. 

Figure  8  shows  a  coaplete  analysis  of  TRW  run  71A  (DREP 
12-3).  The  top  trace  (in  dots)  is  the  current  aagnitude 
along  ship  track  as  aeasured  by  the  ship's  bow  current 
aeter.  In  an  aabient  sea  without  internal  waves,  this  is 
the  ship's  velocity  relative  to  the  water  aass,  which  aay 
itself  be  aoving,  e.g.,  as  in  a  tidal  current.  When  an 
internal  wave  is  present,  its  associated  current  (at  1  a 
below  surface)  will  be  superiaposed  onto  the  ship's  velocity 
relative  to  water,  as  is  apparent  in  this  figure.  The  trace 
close  to  the  current  is  the  Bean  Doppler  shift  aeasured  by 
the  X-band  radar,  converted  into  velocity  units.  It  shows 
the  velocity  of  the  aicrowave  scatterer  relative  to  the 
ship.  The  other  traces,  L-band  Bragg  slope,  X-band  Bragg 
slope,  aicrowave  backseat tered  power,  polarization  ratio, 


video  power,  have  been  explained  above.  The  starting  and 
stopping  of  the  run,  the  polarization  of  the  Microwave  power 
displayed  and  the  incidence  angle  of  the  Microwave  are 
indicated  on  the  top  of  the  figure.  The  direction  of  ahip’e 
heading,  internal  wave  heading  and  wind  vector  are  indicated 
in  the  circled  figure. 

The  following  obeervationa  can  be  Bade  about  Figure  8: 

(i)  The  current  trace  ahows  2  huge  aol iton-1 ike  internal 
wave  "buaps".  The  current  Magnitude  of  the  internal 
wavea  ia  aa  large  aa  70  ca/aec.  Thia  run  ia  an  exaaple 
of  a  very  atrong  internal  wave. 

(ii)  "L-band"  Bragg  alope  ahowa  Modulation  peaka  about  4.3 
and  4.8  dB  (peak  over  background,  POB)  at  t  =  90  aec 
and  t  =  435  aec  whereas  "X-band"  Bragg  alope  ahowa 
slightly  larger  Modulations  of  about  5.5  and  6.7  dB, 
respectively.  This  is  contrary  to  Boat  Modulation 
theories  where  Modulation  of  surface  waves  with  X-band 
wavelength  ia  generally  expected  to  be  a  lot  saaller 
than  that  in  L-band. 

(iii)  Modulation  in  Microwave  backacattered  power  is  a  lot 
higher  at  the  2  Modulation  peaks,  which  are  roughly 
15  dB  and  19  dB  (POB).  Otherwise,  the  rest  of  the 
Microwave  power  trace  ia  alaost  identical  to  that  of 
the  "X-band"  Bragg  slope  trace.  If  the  latter  is 
overlaid  on  top  of  the  forner,  it  can  be  seen  that  even 
subtle  changes  in  the  background  arc  reproduced 
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faithfully.  Thia  suggests  that  except  during  the 
modulation  peaks,  scattering  may  be  predominantly  Bragg 
acattering  so  that  modulation  in  microwave  power  is 
identical  to  modulation  in  Bragg  slope  power,  as 
explained  previously.  During  the  modulation  peaks, 
another  scattering  mechaniam  other  than  Bragg 
scattering  is  responsible  for  the  huge  modulation. 

(iv)  The  polarization  ratio  (hh/vv)  trace  stays  roughly 
constant  at  about  0.25  (or  -6  dB)  except  during  the 
modulation  peaks.  This  is  consistent  with  Bragg 
scattering  which  predicts  a  ratio  of  -6.5  dB.  During 
the  modulation  peaks,  the  polarization  ratio  approaches 
1,  suggesting  specular  reflection  as  the  scattering 
mechanism. 

(v)  The  video  power  trace  confirma  that  specular  reflection 
is  indeed  the  scattering  mechanism  during  the 
modulation  peaks. 

(vi)  One  prominent  feature  in  the  "X-band"  Bragg  slope  is  a 
modulation  null  of  roughly  11  dB  at  t  =  355  sec.  This 
null  is  also  present  in  the  microwave  power  and  to  a 
lesser  extent,  in  the  "L-band"  Bragg  slope.  Since  this 
deep  null  is  detected  by  the  slope  gauge  and  the  radar 
which  are  independent  instruments,  it  is  real  and  not 
an  instrument  malfunction.  Since  we  cannot  correlate 
it  with  any  special  feature  in  the  current,  we  suspect 
its  origin  is  not  solely  hydrodynamic  but  surface  slick 


induced. 


(vii)  The  interpretation  of  the  Doppler  velocity  is 

potentially  the  aost  interesting  but  also  the  aost 
uncertain  in  this  experiaent.  The  reason  for  the 
uncertainty  is  that  there  is  no  other  surface 
instruaent  that  aeasures  surface  scatterer  velocity  to 
confira  what  we  think  the  radar  is  seeing.  In  our 
previous  laboratory  wavetank  experiment* 1 » * >  where  a 
scanning  laser  slope  gauge  is  available  and  where 
correlation  between  successive  scans  is  possible,  we 
found  that  electroaagnetically ,  when  the  radar  is 
seeing  Bragg  scattering,  hydrodynaaically,  there  can  be 
three  different  kinds  of  Bragg  scatterers.  At  low  wind 
speed,  the  Bragg  scatterer  is  a  freely  propagating 
Fourier  component  of  a  slightly  rough  surface  on  top  of 
the  wind  drift  layer,  i.e.,  Doppler  velocity  =  Bragg 
wave  phase  speed  +  wind  drift  layer  velocity.  At 
higher  wind  speed,  the  Bragg  scatterer  could  be  either 
(i)  a  turbulent  patch  aoving  at  the  velocity  of  the 
wind  drift  layer,  i.e.,  Doppler  velocity  -  wind  drift 
layer  velocity  or  (ii)  a  Fourier  component  of  the 
parasitic  capillaries  which  are  phase-locked  to  the 
crest  of  short  gravity  waves,  i.e.,  Doppler  velocity  = 
sh^rt  gravity  wave  phase  speed.  At  high  wind  speed, 
the  Bean  Doppler  velocity  therefore  lies  soaewhere 
between  the  wind  drift  layer  velocity  and  the  phase 
velocity  of  the  doainant  short  gravity  wave.  We  will 
now  consider  these  two  possibilities  in  Georgia  Strait. 


If  the  Bragg  scatterer  ia  a  freely  propagating  Bragg 
wave,  the  offset  between  the  Doppler  velocity  (as 
■easured  by  the  radar)  and  the  current  (as  Measured  by 
the  bow  current  Meter)  should  be  alMost  constant  and 
quite  predictable.  This  is  because  the  offset  would  be 
equal  to  the  Bragg  wave  intrinsic  phase  speed  (which  ia 
constant)  +  the  wind  drift  layer  velocity  (which  is 
usually  relatively  smsII  and  alMost  constant  if  the 
wind  field  is  alaost  constant).  If  the  Bragg  scatterer 
is  a  combination  of  turbulent  patches  and  parasitic 
capillaries  of  short  gravity  waves,  then  the  offset 
between  the  Doppler  velocity  and  the  current  would  be 
More  variable.  This  is  because  the  offset  would  be 
equal  to  the  aean  speed  of  turbulent  patches  and  short 
gravity  waves  in  the  field  of  view.  The  speed  of  the 
turbulent  patches  would  be  the  wind  drift  layer  speed 
which  is  quite  constant  and  predictable.  The  speed  of 
short  gravity  waves  (More  precisely,  its  coMponent 
along  ship  track)  would  be  More  variable  due  to 
variability  in  both  Magnitude  and  direction.  Referring 
now  to  Figure  8,  we  can  see  that  the  offset  between 
Doppler  velocity  and  current  Meter  output  is  far  from 
constant.  This  rules  out  that  scattering  is  entirely 
due  to  Bragg  scattering  froM  freely  propagating  Bragg 
waves.  Closer  exaaination  will  show  that  the  Microwave 
scatterers  are  Moving  towards  the  ship  before  and 
between  the  internal  wave  "buMps”,  but  during  the 


buaps"  they  actually  change  direction  and  love  in  the 
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•aie  direction  as  the  ship  and  internal  wave.  To 
understand  how  this  happens,  we  refer  to  Figure  9 
which  characterizes  alaost  all  the  natural  internal 
waves  we  have  seen  in  Georgia  Strait.  Based  on  visual 
observations,  an  internal  wave  propagating  froa  left  to 
right  is  usually  preceded  by  a  dead  cala  region.  Then 
coaes  a  rough  region.  The  deaarkation  between  the  cala 
and  rough  region  is  extreaely  sharp.  Froa  the 
deaarkation  line,  short  gravity  waves  appear  to  be 
growing  in  both  wavelength  and  aaplitude  as  they  travel 
to  the  right,  in  the  saae  direction  as  the  internal 
wave,  irrespective  of  wind  direction.  The  short 
gravity  waves  aay  or  aay  not  break  at  soae  point, 
depending  on  whether  the  internal  wave  is  strong 
enough.  Further  to  the  right,  the  short  gravity  waves 
grow  longer  but  slowly  die  down  in  aag-nitude,  becoaing 
another  dead  cala  region  in  front  of  the  next  cycle. 

In  Run  71A  as  depicted  in  Figure  8,  the  wind  is 
blowing  at  ~  115°  whereas  the  internal  wave  is 
traveling  at  354°.  Outside  the  internal  wave  buaps, 
the  aicrowave  scatterers,  irrespective  of  origin,  are 
naturally  aligned  with  the  wind  direction,  so  that  the 
net  coaponent  along  ship  track  is  towards  the  ship. 
Inside  the  internal  wave  buaps,  the  short  gravity  waves 
as  depicted  in  Figure  9  are  now  the  aicrowave 
scatterers  and  so  travel  in  the  saae  direction  as  the 
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internal  waves  and  ship.  In  factr.  some  of  the  sharp 
excursions  downwards  in  the  Doppler  velocity  can  be 
clearly  identified  with  the  peaks  in  microwave  power 
and  peaka  in  polarization  ratio.  This  closely  agree* 
with  our  notion  that  specular  reflection  coses  f roa 
short  gravity  waves  breaking  and  surging  forward. 
Outside  the  internal  wave  buaps,  we  have  shown  that 
polarization  ratio  indicates  the  scattering  to  be 
likely  Bragg  scattering.  The  Doppler  velocity  itself 
however  cannot  tell  us  whether  the  Bragg  scatterers  are 
(i)  freely  propagating  or  (ii)  a  combination  of 
turbulent  patches  and  parasitic  capillaries  on  short 
gravity  waves.  If  (i)  is  correct,  the  offset  between 
the  Doppler  velocity  and  current  traces  should  be  equal 
to  Bragg  wave  intrinsic  phase  speed  (24  ci/sec)  +  wind 
drift  velocity  along  ship  track  (3.8  s/sec  x  0.03 
x  cos  66°  =5.1  ci/sec)  =  29.1  ci/sec.  This  is  not  in 
bad  agreement  with  the  figure.  However,  (ii)  cannot  be 
ruled  out.  In  fact,  a  combination  of  (i)  and  (ii)  can 
also  occur  and  cannot  be  ruled  out.  Measuring  the 
offset  between  Doppler  velocity  and  current  can  only 
definitively  show  whether  (i)  is  incorrect,  e.g.,  if 
the  offset  is  60  ci/sec,  then  it  cannot  possibly  be  due 
to  a  freely  propagating  Bragg  wave  on  top  of  a  wind 
drift  layer  in  this  low  wind.  A> scanning  laser  slope 
gauge  is  definitely  needed  to  shed  sore  light  on  the 
identity  of  the  Bragg  scatterers. 


(viii)  In  Run  71A  in  Figure  8,  the  debarkation  between  the 
calb  region  and  rough  region  in  the  internal  wave 
occura  at  the  divergence  zone  in  the  current.  In  all 
the  other  runs  we  have  analyzed,  thia  debarkation 
happens  at  the  convergence  zone.  Juat  in  caae  there  is 
any  doubt,  we  have  analyzed  Run  69  which  is  the  ship 
going  south  one-half  hour  earlier  and  crossing  the  sane 
internal  waves.  The  debarkation  between  rough  and  calb 
regions  still  lie  in  the  divergence  zone.  This 
peculiarity  is  still  a  systery  to  us. 

(ix)  The  "X-band"  Bragg  slope  and  "L-band"  Bragg  slope  have 
been  analyzed  by  tracking  the  "Bragg  frequency"  with 
the  Doppler  velocity  of  the  radar.  We  have  found  that 
tracking  with  the  current  or  a  constant  ship  velocity 
does  not  bake  buch  difference.  This  ibplies  that  the 
effect  of  bodulation  is  a  lot  greater  than  whatever 
error  we  bay  have  incurred  by  using  a  slightly  off 
Bragg  frequency. 

We  will  now  exabine  TRW  Run  81B  which  is  an  exabple  of  a 
very  regular  internal  wave  train  (Figure  10,  see  Caponi  et 
al,  this  issue,  for  a  description  of  the  conditions  of  TRW 
Run  81B).  We  can  bake  the  following  observations: 

(i)  The  agreebent  between  the  "X-band"  Bragg  slope  and  the 
bicrowave  power  is  even  bore  ibpressive  in  this 
exabple.  When  the  slope  trace  is  overlaid  on  top  of 
the  power  trace,  it  can  be  seen  that  the  only 
disagreeaent  occurs  around  the  null,  where  the  noise 
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level  of  the  elope  gauge  has  been  reached  (and  thus 
shows  apparently  higher  level)  and  around  the 
■odulation  peaks,  where  spikes  froa  specular  reflection 
occur  in  aicrowave  power  but  not  slope. 

(ii)  Coaparison  between  the  "L-band"  slope  and  the  X-band 
■icrowave  power  will  show  that  while  there  is 
siailarity,  the  agreeaent  is  nowhere  as  close.  If  the 
whole  slope  spectrua  is  integrated  to  give  the  aean 
square  slope,  it  can  be  shown  that  the  agreeaent 
between  that  and  the  aicrowave  power  is  even  worse. 

This  shows  that  "X-band"  Bragg  slope,  as  we  have 


defined  it,  is  the  proper  descriptor  of  the  surface  to 
be  coapared  with  X-band  aicrowave  backscattering  power, 

(iii)  The  direction  of  travel  between  internal  wave  and  the 
ship  can  be  inferred  froa  this  figure.  Referring  back 
to  Figure  9,  we  can  see  that  the  ship  is  first 
encountering  the  rough  region,  then  the  very  sharp  drop 
to  the  cala  region.  This  aeans  the  ship  is  travelling 


against  the  internal  wave. 

(iv)  The  polarization  ratio  is  again  close  to  0.25  (-  6dB) 


aost  of  the  tiae,  indicating  the  presence  of  Bragg 
scattering.  However,  other  than  peaks  that  correspond 
to  real  power  peaks,  there  are  also  spurious  peaks  in 
the  aodulation  nulls.  This  is  an  artifact  because  both 
vv  and  hh  signals  are  dropping  close  to  noise  level. 


i  livi.  r»3*».r 


(v)  The  sharp  debarkation  lines  between  rough  and  cals 

regions  occur  at  the  convergence  zone  of  the  current  in 
this  run. 

We  will  now  examine  TRW  Run  112  (Figure  11).  This  is  the 
case  of  a  very  weak  internal  wave  in  a  call  sea.  The 
following  notable  features  are  particular  to  this  run: 

(i)  Even  though  the  internal  wave  is  weak,  aodulationa  in 
slope  and  aicrowave  are  nonetheless  quite  large, 
exceeding  20  dB  in  aicrowave.  Again  the  noise  floor  in 
the  slope  masks  the  real  bottoas  of  the  aodulation 
cycles . 

(ii)  If  we  examine  carefully,  almost  every  convergence 
point  in  the  current  corresponds  to  the  sharp  dividing 
line  between  the  cala  region  and  rough  region  of  an 
internal  wave  cycle.  The  aost  prominent  ones  are 
pointed  out  by  arrows  in  the  figure. 

(iii)  The  wind  is  extreaely  low  (1.9  a/s),  the  sea  very 
cala  and  the  internal  wave  weak  in  this  run.  This 
would  seea  like  a  good  candidate  for  observing  freely 
propagating  Bragg  waves.  Indeed,  the  Doppler  velocity 
does  show  an  alaost  constant  offset  of  roughly 

35  ca/sec  froa  the  current.  A  freely  propagating  Bragg 
v««e  would  have  a  constant  offset  of  roughly  28  ca/sec. 
The  conjecture  therefore  seeas  quite  reasonable. 


5. 


CONCLUSIONS 


With  respect  to  specular  reflection,  we  can  draw  the 

following  conclusions: 

(i)  Specular  reflection  is  highly  significant  at  20° 
incidence  angle  and  not  negligible  at  40°. 

(ii)  Our  relatively  crude  instrumentation,  which  was  meant 
at  first  for  detection  purposes  only,  indicates  that 
the  specular  points  have  large  enough  radius  that  they 
are  no  longer  Rayleigh  scatterers  at  X-band  and 
therefore  that  the  specular  points  observed  by  the 
optical  system  are  also  specular  to  the  X-band 
microwave  radar. 

With  regard  to  comparison  of  microwave  scattering  with 

slope  data,  we  can  draw  the  following  conclusions: 

(i)  The  good  agreement  of  the  X-band  Bragg  slope  and  the 
microwave  backscattered  power  (when  specular  reflection 
is  absent)  shows  that  both  the  DRBP  slope  gauge  and  the 
TRW  radar  are  highly  accurate  instruments  over  their 
respective  dynamic  ranges.  It  also  shows  that  Bragg 
scattering  is  a  good  description  of  microwave 
scattering  when  specular  reflection  is  absent. 

(ii)  X-  and  L-band  modulations  (as  evidenced  by  the 
respective  Bragg  slopes)  are  comparable  in  most  runs. 
This  is  quite  unexpected  for  most  modulation  theories. 

(iii)  In  calm  seas,  lower  signal  level,  but  higher 
modulations  are  observed.  These  modulations  are 
primarily  due  to  Bragg  scattering.  In  rougher  seas, 


the  hydrodynamic  aodulationa  are  generally  lower.  The 
microwave  aodulationa  car.  however  be  significantly 
enhanced  by  specular  reflection. 

(iv)  Deep  minima  have  been  observed  near  current  patterns, 
but  with  no  obvious  relationship  to  the  current 
variations.  These  may  be  due  to  the  presence  of 
organic  films  on  the  surface. 

(v)  Under  calm  seas,  weak  internal  waves  can  produce  very 
strong  modulations  at  X-band  but  weaker  modulations  at 
L-band  (as  evidenced  by  their  respective  Bragg  slopes). 
Again,  organic  film  is  a  likely  non-hdyrodynamic  factor 
that  should  be  checked  into  in  future  experiments. 

Our  previous  laboratory  work<1,2>  suggests  that  in 
addition  to  Bragg  scattering,  wedge  diffraction  from  short 
gravity  wave  crests  and  specular  reflection  from  breaking 
waves  are  important  mechanisms  for  microwave  backscattering 
at  moderate  incidence  angles  (20°-70°).  Generally  speaking, 
wedge  diffraction  and  specular  reflection  enhance  scattering 
by  0(3-4  dB)  and  0(10  dB)  respectively.  To  discern  the 
wedge  diffraction  component  in  ocean  scattering,  a  scanning 
laser  slope  gauge  needs  to  be  co-located  with  the  radar  so 
that  a  deterministic  study  can  be  carried  out.  In  the 
present  experiment,  wedge  diffraction  nay  have  enhanced  the 
Bragg  acsttering  background  quite  uniformly  by  a  few  dB.  It 
may  also  be  buried  around  or  under  the  specular  reflection 
events.  (Wave  crests  tend  to  become  sharp  wedges  Just 


before  they  break.)  The  confirmation  of  wedge  scattering  in 
the  ocean  will  have  to  wait  for  a  more  precise  measurement. 
The  video  instrumentation  in  this  experiment  were  meant 


A 


to  detect  specular  facets  and  provide  a  rough  estimate  of 
their  sizes.  A  quantitative  and  accurate  measurement  of  the 
radius  of  curvature  distribution  of  specular  points  on  the 
ocean  surface  will  require  improved  instrumentation. 
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APPBNDIX 


A  specular  point  in  the  Microwave  antenna  footprint 
will  appear  as  a  white  dot  or  patch  in  the  video  recorded 
image.  This  is  due  to  the  fact  that  the  finite-size 
diffused  laap  has  produced  a  finite-size  virtual  image  in 
the  convex  surface  around  the  specular  point.  From  simple 
geometric  optics,  the  magnification  is  given  by 

m  *  r  +  2*0  ( 1 ) 

where  r  is  the  radius  of  curvature  at  the  specular  point  and 
0  is  the  object  distance  (i.e.,  laap  to  water  surface 
distance) .  The  real  length  of  the  virtual  image  is  given  by 

ti  =  *o.m  (2) 

where  i-i  =  length  of  the  image  and  S-o  =  length  of  the 
object,  i.e.,  O.D.  of  the  laap. 


then 


The  apparent  length  of  the  image  in  the  video  camera  is 


0  % 

0  m  Q  •  -  —  —  m  _ 

ia  i  0  +  i  1  +  m 


where  ha  -  apparent  length  of  the  image  in  the  video  camera 
and  i  =  distance  of  virtual  image  behind  the  water  surface. 

Combining  Bquations  (1),  (2)  and  (3),  the  radius  of 
curvature  can  be  expressed  in  terms  of  the  apparent  length 
of  the  image  in  the  video  camera  as: 


4 

te)- 


We  note  here  that  if  a  projection  laap  is  used  instead  of  a 
diffused  laap,  the  filament  of  the  laap  will  become  the 
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object  instead,  the  virtual  iaage  will  be  a  point  source, 
and  no  aize  differentiation  will  be  possible.  The  intensity 
of  the  virtual  iaafe  will  be  relatable  to  the  radius,  but 
the  intensity  seasuresent  of  a  point  is  obviously  less 
accurate  than  the  size  aeasureaent  of  a  patch  in  a  video 


Figure  5 


Size  distribution  of  specular  facets  In  TRW  Run  81B  as  a  function 
of  kr  In  both  X-  and  L-bands.  20  sec.  of  data;  40°  Incidence 
angle;  2.8  m/s  wind  speed. 
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Figure  6.  Power  distribution  of  specular  facets  in  TRW  Run  81B  as  a  function 
of  kr  in  both  X-  and  L-bands.  20  sec.  of  data;  40°  incidence 
angle;  2.8  m/s  wind  speed. 
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